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ABSTRACT: There is a newly emerging understanding that in the chaotic domain of isolated 
finite interacting many particle systems smoothed densities define the statistical description of these 
systems and these densities follow from embedded (two-body) random matrix ensembles and their 
various deformations. These ensembles predict that the smoothed form of matrix elements of a 
transition operator between the chaotic eigenstates weighted by the densities at the two ends (i.e. 
the bivariate strength density) will be a bivariate Gaussian with the bivariate correlation coefficient 
arising out of the non-commutability of the hamiltonian and the transition operator involved. The 
ensemble theory extends to systems with a mean-field and a chaos generating two-body interaction 
(as in nuclei, atoms and diffusive quantum dots). These developments in many-body quantum 
chaos are described with special reference to one-body transition operators. 

1. Introduction 

Random matrix physics, started with Wigner and Dyson's GOE, GUE and GSE hamiltonian ran- 
dom matrix ensembles Q , has grown far and wide in the last fifteen years with the investigation of 
a variety of random matrix ensembles (interpolating, banded, partitioned, chiral etc.) and applica- 
tions in almost all branches of physics Q (most recent application being in Econophysics Q). As 
Guhr et al Q state 'the universal validity of random matrix theories is reminiscent of the success 
of thermodynamics in the last century • • •'. Recently it is recognized by large number of research 
groups in atomic, molecular, nuclear and mesoscopic physics (see for example ||, [| |§> 0])> m the 
context of many-body quantum chaos, that embedded random matrix ensembles (EE) (in particular 
EGOE(/c), the embedded Gaussian orthogonal ensemble of random matrices of fc-body interactions) 
are relavant for 'isolated finite interacting particle systems'. Our purpose in this article is to discuss 
sum of the results of EE for transition matrix elements. In Section 2 definition and construction 
of EGOE(2) for fermion systems are given; introduced here are also various deformed EGOE. The 
EGOE(k) results for density of states and transition matrix elements are briefly described in Sec- 
tion 3. The topic 'quantum chaos and transition strength sums' form Section 4. In Section 5, 
EE(l+2) results for matrix elements of one-body transition operators with emphasis on the role of 
the mean-field basis and strength functions are presented. Section 6 gives concluding remarks. 



2. Embedded Ensembles 



EGOE(/c) for many fermion systems is defined by assuming that the many particle space is a 
direct product space, of single particle states, as in the nuclear/atomic shell model. EGOE(fc) 

invited talk in the National Conference on 'Dynamical Systems: Recent Developments ' held at Unversity 
of Hyderabad, Hyderabad during 4-6th November, 1999. 
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for m (m > k) fermion systems (with the particles distributed say in N single particle states) is 
generated by defining the hamiltonian H, which is say /c-body, to be GOE in /c-particle space and 
then propagating it to m - particle spaces by using the geometry of the m - particle spaces H . Let 
us consider EGOE(2). Given the single particle states \v%) , i = 1,2, ... ,N, operator form of the 
2-body hamiltonian is (with a in ( | | ) indicating that the states involved are antisymmetrized) , 

H= Yl {vkVi\H\vi Vj ) a a\ at k a n a Vj (1) 

Symmetries for the two-body matrix elements (TBME) (u k v\ \ H \ i/j Vj) in (1) are, 

(u k vi\H\ Vj Vi) a = - (u k ui\H\ui Vj) a , (u k ui\H\ui Vj) a = (u { Vj \ H \ v k u{) a (2) 

The hamiltonian H in m-particle spaces is defined in terms of the TBME via the direct product 
structure. The non-zero matrix elements of the m-particle H matrix are of three types, 

(VXV 2 ■■■V m \H\ V X Vl ■ ■ ■ V m ) a = Y I H | ViVj) a 

Vi<Vj<V m 



(v p V 2 Vz ■ ■ ■ U m | H I V\V % ■ ■ ■ V m ) a = Y Wp"i I H I V\Vi) a 
(u p V q U 3 ■ • • V m I H | VlV 2 V 3 ■ ■ ■ V m ) a = (VpV q \ H I V\V 2 ) a 

EGOE(2) is defined by (l)-(3) with GOE representation for H in two-particle space, 
(v k v\ | H | Vi Uj) a are independent Gaussian random variables 



(3) 



(4) 



(u k vi | H | i/j Uj) a = 0, \(v k n | H | z/j Vj) a \ = v\l + %),(fei)) 



In (4) bar denotes ensemble average and v is a constant. Note that the H(m) matrix dimension d is 
d(N,m) = ( ^ ) an d the number of independent matrix elements (ime) are ime(N) = d 2 (d 2 + l)/2 
where the two-particle space dimension d 2 = N(N — l)/2. For example, <i(12,6) = 924 and 
ime(12) = 2211. The EGOE(2) is also called TBRE (two-body random matrix ensemble). Exten- 
sion of (1-4) for boson systems is straightforward and it is described elsewhere [^]. Hamiltonians 
for many interacting particle systems contain a mean-field part (one-body part h) and a two-body 
residual interaction V mixing the configurations built out of the distribution of particles in the 
mean-field single particle orbits; h is defined by single particle energies (SPE) ej, i = 1 — N 
and V is defined by TBME. Then it is more realistic to use EE(l+2), the embedded ensemble of 
(l+2)-body hamiltonians, 

EE(l+2) : {H} = [h(l)]+X{V(2)} . (5) 

Here {V} is EE(2), i.e. it is EGOE(2) with v = 1 in (4) or an ensemble with TBME being 
independent random variables with a distribution different from Gauassian (for example uniform 
distribution). Similarly [h] is a fixed hamiltonian or an ensemble with SPE chosen random but 
following some distribution. Finally, [h] and {V} are independent. It is to be expected that the 
generic features of EE(l+2) approach those of EGOE(/c) for sufficiently large values of A and 
significant results emerge as A is varied starting from A = 0. Let us mention that EE(l+2) is also 
called TBRIM (two-body random interaction model). A second class of EE are the partitioned 
embedded ensembles (p-EE) where the hamiltonian in two-particle space is block structured (i.e. 
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the space divides into subspaces ©r) with variances of the matrix elements in each block being 
v^ T , for the block connecting T and T' subspaces. Third class of EE are EE-sym where v^/ = 
for r / T' and these ensembles are important for Hamiltonians carrying symmetries (for example 
J or JT in nuclei and J or LST in atoms). Finally there are also the modified X+EE ensembles 
where K is a fixed operator. 



3. Basic results for EGOE(k) 



Generic results for state densities and transition matrix elements, for EGOE(k), that are essentially 
valid in the dilute limit, which correspond to N,m,k — > oo, m/N — > and k/m — > 0, are well known. 
The eigenvalue density 1(E) or its normalized version p(E) takes Gaussian form |5|, 

1(E) = ((S(H -E)))=d p(E) ; p(E) E °°? pjE) = P g(E) = ^exp - \ {^) ' (6) 

The binary correlation approximation, originally used by Wigner for deriving the semi-circle state 
density for GOE (k = m in EGOE(k)) is used by Mon and French Q to derive (6) via the m- 
particle space moments (H v ) m of 1(E). Firstly it is seen that by definition all odd moments of 
1(E) will vanish. Using the normalization that (H 2 ) m ~ 2 = 1 (then (H 2 ) m = ^ ™ J) one has the 

basic result H(k)o(t)H(k) => ^ m ~ * ^ o(t) under ensemble average in m-particle spaces and in the 

dilute limit. Using this and that in the trace {H p ) m binary associations dominate, one can drive 
formulas for the moments by writing down all possible binary association diagrams. For example 
for (H 4 ) there are three diagrams, AABB ABBA © ABAB => 2{AABB} © ABAB. Note that 
A and B are -ff-operators and they are fc-body in nature. Evaluating the irreduciable diagrams 
AABB and ABBA give the 4th reduced moment ^4 and the 4th cumulant k^, 

Similar results for the 6th and 8th cumulants are derived: ke = k 3 (6k — l)/m 2 + 0(l/m 3 ) and 
k 8 = -Ak 5 (23k-9)/m 3 +0(l/m 4 ). Thus in the dilute limit one recovers the Gaussian form for state 
densities (note that we need in fact not k/m — > but k 2 /m — > 0). Thus, for a two-body interaction 
m ~ 12 gives a good Gaussian. Note that for m = 4 one has k/t = — 1 implying semi-circle shape as 
seen in many numerical calculations. In practice one has to apply Edgeworth (mostly 3rd and 4th 
moment /cumulant) corrections to the Gaussian form. More important point, though not discussed 
here, is that local level fluctuations given by EGOE(k) are of GOE type. 

Similar to (6), for EGOE(k) the transition matrix elements weighted by the densities at the two 
ends (i.e. bivariate strength densities) take bivariate Gaussian form flOtl , 



I bi v,o(E, E>) = ((C*6(H - E')05(H - E) f 
I'(E')\(E' I O I E)\ 2 I(E) = ((&0)) p bi v,o(E,E') 



Pbiv,o(E, E') E 92^ p biv . (E, E') — Pbiv-g-, o(E, E' ; ei, €2, &i, &2, C) 



1 



exp ■ 



2na 1 a 2 V^C 1 I 2 ( l ~ C 



E-e!\ 2 / g-eA / E'-e 2 \ | (E'-e^~ 



&1 J \ V\ J \ &1 J \ 02 



(8) 
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The bivariate reduced central moments of puv;0 are \i vq = (o^ ( H a ^' 2 \ Q O ( ^Ff 1 ) ) / {p^O\ and 
C = /in is the bivariate correlation coefficient. Let us consider the evaluation of fi pq by representing 
H by EGOE(k) and the transition operator O by EGOE(t). It is also assumed that H and O 
ensembles are independent. Now the correlations in [x vq arise due to the non-commutability of H 
and O operators. Firsly it is seen that all fi pq with p + q odd will vanish on ensemble average and 

also fi pq = n qp . Moreover a\ = a\ = (o^OH 2 ^j / l^O^O^j = ( 'fc )• The ^ rs ^ non-trivial moment 
Mil is > 

C = Mil = (0^{t)H{k)0{t)H{k)) m I Uoi(t)0(t)\ m (H{k)H(k)) m } -► ABAB/(i4^)(BB) 

fct fc(ib - l)t(t - 1) 



m — t 
k 



-i 



m 
k 



1 



777 



+ 



2m? 



+ 0{l/m i 



(9) 

Let us now consider the cases with p + q = 4. The diagrams for these follow by putting and O 
at appropriate places in the (H 4 ) diagrams, 



= 2 

Ml3 



m — k \ ( m 
k )\k 

(O^AABOB) © (O^ABBOA) © (0+ AB AOB)} / [(CfiO) {(H 2 )) 2 } 
m — t \ f m \ fm — k\fm — t \ ( m 

k ) V k ) V fc / V fc /J v fc 

M22 = [(O^AAOBB) © (O^ ABOBA) © (OUSOAB)] / [(0t0) ((tf 2 )) 2 
m — t \ I m — k — t \ I m — t 



k 0i 
h3 



k J ' \ k 
k 40 = /i4o — 3 = 

hi = M31 - 3/in = C&04 



m — k 
k 



rn 
k 



777, 



777 
k 



fc 2 (fc-l) 2 

2m 2 



+ 0(l/m 3 ) 



(10) 



k 22 = M22 - 2C 2 - 1 = C 

k 2 k 2 [(k - l) 2 + Akt - 2t] 



m — k — t 
k 



m — t 
k 



in 



2m 2 



+ 0{l/m 3 ) 



Further results and numerical tests etc. are given in @. Finally it is worth mentioning that the 
EGOE smoothed Gaussian forms for the state densities (6) and the bivariate strength densities (8) 
together with the result that local strength fluctuations for EGOE (just as for GOE) are of Porter- 
Thomas type, are used recently to derive the EGOE formulas for information entropy (S in f°) and 
number of principal components (NPC), which are measures of complexity and chaos, in transition 
strength distributions 

4. EE(l-)-2) and quantum chaos: transition strength sums 



The principal reason for EGOE receiving considerable attention in the last 3-4 years is due to the 
fact that: there is a newly emerging understanding that in the chaotic domain of isolated finite 
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interacting many particle systems smoothed densities (they include strength functions) define the 
statistical description of these systems and these densities follow from embedded random matrix 
ensembles. This conjucture is based on the calculations for transition strength sums 
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Fig. 1. (a) Occupation numbers for a 25 member EE(l+2) ensemble, defined by the hamiltonian 
h(l) + X{V(2)}, in the 924 dimensional N = 12, m = 6 space. Results are shown for the lowest 5 
single particle states and for six values of A. In the calculations occupation numbers are averaged 
over a bin size of 0.1 in (E — e)/o~; e is centroid and a is width. The spectra of all the ensemble 
members are first zero centered and scaled to unit width and then the ensemble average is carried 
out. The estimate of [11] gives A c ~ 0.05 for order-chaos border in the present EE(l+2) example. 



It is clearly seen that once chaos sets in, the occupation numbers take stable smoothed forms. For 
A = .08, the smooth forms are well represented by ratio of Gaussians (smooth curves in the figure), 
(b) same as (a) but for the transition strength sum generated by the one-body transition operator 
a\a$. The bin size used here is 0.01. 



carried out using different types of interpolating hamiltonians generating order-chaos transitions 
in many different systems: (i) for occupancies using a 20 member EE(l+2) in 330 dimensional 
N = 11, m = 4 space with h(l) defined by the single particle enegies e$ = £ + (1/i) ; i = 1, 2, . . . , 11 
and V(2) is EGOE(2) |l2[]; (ii) for occupancies using the four interacting electrons Ce atom Q; 
(iii) for occupancies using a symmetrized coupled two-rotor model |L3|]; (iv) using two different 
types of interpolating ensembles in nuclear shell model for 24 Mg for the Gamow- Teller strength 
sums and for occupancies ||. Most significant conclusion of all these studies is that transition 
strength sums (note that occupancies are one-particle transfer strength sums) show quite different 
behaviour in regular and chaotic domains of the spectrum. In order to make this argument clear, 
calculations are carried out using EE(l+2) as in (i) but for iV = 12 and m = 6. The EGOE gives 
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via (8) that the transition strength sum density I ^ (E) = ((O^OS(H — E)jj, which is a marginal 

density of the bivariate strength density, is a Gaussian. Therefore, using (6), it is immediately seen 
that transition strength sums vary with excitation energy as ratio of Gaussians. Fig. la shows 
results for occupancies and Fig. lb for transition strength sums for the one-body operator a} 2 a§ 
calculated for various values of the interpolating parameter A in the EE(l+2) hamiltonian. From 
Figs. la,b it is clearly seen that below the region of onset of chaos transition strength sums show 
strong fluctuations (in the regular ground state domain perturbation theory applies). In this region 
there is no equilibrium distribution for Iq^q(E), 1(E) and other densities. However in the chaotic 
domain (there are methods for determining the critical A c that marks order-chaos border 0, |Tl| ) 
the densities can be replaced by their smoothed forms. Note that fluctuations in strength sums are 
basically given by 1/NPC. Therefore there is a statistical mechanics, defined by various smoothed 
densities and they are given by EGOE, operating in the quantum chaotic domain of isolated finite 
interacting particle systems (this is also the essence of statistical nuclear spectroscopy ju], [l4|, 15] 



and also statistical spectroscopy in atoms [||]); see also Q. In fact in favourable situations, it is 
possible to introduce thermodynamic concepts (effective temperatures and chemical potentials etc.) 
in the chaotic domain [Ej|. 

5. EE(l-)-2) results for matrix elements of one-body transition operators: role of the 

mean-field basis and strength functions 

5.1 Strength functions: transition from Breit - Wigner to Gaussian form 

Decomposing the m-particle spaces into subspaces T (say defined by the irreps of a group structure 
in the Hilbert space defined by (N,m)), i.e. m — » I\ gives, 

1(E) = ((5(H = £ ((5(H - E))f = £ f(E) ; f (E) — > T(E) = I T g{E) 



I K (E) = ((K5(H - E))) = £ ((K5(H - E))f = ; I T K (E) — I?(E) = f K , g (E) 

r r 



(11) 

There are plausible arguments in favour of the Gaussian forms in (11) for EGOE. The decompo- 
sitions in (11) are related to the well known 'strength functions' (also called local spectral density 
of states (LDOS) in literature) which are basic ingredients of a many particle system. Given a 
compound state (j>k, the probability of its decay into stationary states i/je (generated by H) is given 
hy \{4>k I V'b)) 2 - Then the strength function F^(E) is, 

\4>k)=Y, C k \1>e) ; F k (E)=J2\ci'\ 2 S(E-E') = {5(H-E)) k (12) 

E E> 

If we consider a set of states (f>k that belong to a irrep T as in (11), then the average strength 
function Fr(E) is nothing but the partial densities Ir(E) in (11); Ir(E) = J2ker Ek(E) and Fr(E) = 
Ir(E)/d(T). A quite different and useful way to look at strength functions is to think of 4>k as a 
compound state generated by the action of a transition operator O on a state ipE k (for example 
ground state), 

\<PE k ) = T7> ^ \<PE k <PE k ) = 1 

((0U(H - E)05(H - E k )W ( U,) 

\M = E cl |^) * f, M e) = " {{omH _ Ek))) ' 
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Fig. 2. Strength functions (Fs k (E)) for the EG0E(l+2) ensemble defined in Fig. 1. The \k) 
states are the mean-field h(l) states defined by the distribution of m particles in the ./V single 
particle states and their energies are = (k\h(l) + XV(2)\k). In the calculations E and are 
zero centered for each member and scaled by the spectrum (E's) width a; E^ = — e)/a and 
E = (E — e)/<7. For each member |C^| 2 are summed over the basis states \k) in the energy window 
Eft db A and then ensemble averaged Fg, (E) vs E is constructed as a histogram; the value of A is 
chosen to be 0.025 for A < 0.1 and beyond this A = 0.1. Results are shown for A = 0.01, 0.05, 0.1, 
0.2, 0.3 and 0.5. The histograms are EGOE(l+2) results, continuous curves are BW fit and dotted 
curves are Edgeworth corrected Gaussians (ED); the ED incorporates |T(| skewness (71) and excess 
(72) corrections. For A = 0.01, 0.05 and 0.1 only BW, for 0.2 both BW and ED and for 0.3 and 0.5 
only ED are shown. In the figure all the results are for E/~ = 0. 



Starting with the EE(l+2), it is to be expected that for sufficiently large values of A in (5), EGOE(fc) 
description should be valid and therefore applying (8) gives the shape of the strength function to 
be Gaussian (conditional density of a bivariate Gaussian is a Gaussian) and hence its width a k is 
independent of the energy E^ or k. These results are indeed seen in numerical calculations. The 
important question is how F^{E) changes as A is varied. The standard form, normally employed 
in many applications, for strength functions is the Breit-Wigner (BW) form characterized by a 
spreading width T^, 

1 r 

F k ; BW (E) = - {E _ Wk) k 2 + Tl/4 (14) 

where E^ = (k | H | k) = § 00 00 Fk(E)E dE. It is expected that BW form and the Gaussian form 
should appear as A is varied and this question is investigated in: (i) using EE(l+2) and constructing 
Fk(E) for various values of A, here <pk are the mean-field (h) basis states (only final result of this 



7 



study is reported in p2]); (ii) by carrying out nuclear shell model calculations with H\ = h+ A V for 
28 Si with (f)^ chosen to be mean- field basis states |l7|| ; (iii) using the three-orbital Lipkin-Meshkov- 
Glick model [jl8| 1; (iv) using a symmetrized coupled two-rotor model Results from these studies 
show that when the system is chaotic (so that level and strength fluctuations follow from GOE), 
strength functions take Gaussian form and order-chaos transition implies for Fk(E) transition from 
BW to Gaussian, 

order — > chaos Fk-.Bw(F) — > Fk-.g(E) (15) 

For further understanding of the nature of Fk(E) we performed EGOE(l+2) calculations in the 
924 dimensional N = 12 and m = 6 space and the results are shown in Fig. 2. From the figure it is 
clearly seen that there is BW to Gaussian transition and the value of the interpolating parameter 
A = \p k for onset of this transition is \p k ~ 0.2. This should be compared with A = A c ~ 0.08 
for onset of chaos in level fluctuations, occupancies and strength sums of one-body operators in 
the present EGOE(l+2) example (Fig. 1). A similar BW to Gaussian transition is seen in [jlT], 
and in the nuclear shell model example considered in |17| A c ~ 0.3 and \F k ~ 0.6. Thus the BW 



form for Fk(E) extends into the chaotic domain (A > A c ) and the transition to Gaussian shape 
takes place in the second chaotic layer defined by A > \F k ■ Full understanding of A p k is at present 
lacking. 

5.2 matrix elements of one-body transition operators in the chaotic domain 

Let us consider a one-body transition operator written in occupation number representation O = 
Ca/3 a a a >/3 an d its matrix elements | (Ef \ O \ Ei)\ 2 in many-particle eigenstates expanded in the mean- 
field basis states \kt) are (with e a being single particle energies and £i being the mean-field basis 
states energies), 



(E f \0\E i )\ 2 = l £cf {k f \0\h 



diag + offdiag 



i Kf 



c. 



\(k f | O I h)\ 2 + E t C H Ck> Cf; C% (k f I O I h) (k' f I O | k\ 



(16) 

The diag term in (16) involves only the strength functions (see (12)) and also for this term, for a 
given kf and k{ only one e a p in O will contribute. Following |15| it is easily seen that, 



m\o\E i )\l iag = Y J K?\ 2 EMi 

a (3 { £, 



n a )Y 



c 



r< E f 



(17) 



Assuming that (np(l — n a )) £l do not vary much over the number of principal components (say 
N e ff is NPC for the mean-field basis states), it can be replaced by its corresponding average (this 
is verified using EE (1+2) in p9| ). Then, 



Ef I O | Et)\* = E l^l 2 (n fi (l - n a )f* E \ C l 



a p 



a 



£f—£i—£f3 + £a 



(18) 



The |C| 2 's in (18) are nothing but strength functions. With the assumption that strength functions 
are a function of (E — £i)/si where Si is a scale parameter (spreading width Yi for BW and the 



spectral width o~i for Gaussian), we have the results (with D{E) denoting mean spacing), 



F k {E) = (sk^fdE - £ k )/s k ) , jF k (E)dE 

D(E)F k (E) ^(NejffiO))- 1 f((E-£ k )/s k 

2 d£ k 



1 



D(E) 



D(E) 
F k (E)d£ k 



WefffiO)} 



-l - 

s 



(19) 



Substituting F's for |C| 2 in (18) and replacing the £{ summation by an integral, both using (19), 
give the final result in terms of occupancies and the mean spacings, 



n, 



\(E f \0\E i )\l iag = Y,\e«p? (Ml 

= E M 2 M 1 " ^)>^ D{E f ) T(A, s h s f ) ; 

a,/3 

i iy+iy 



)) Si £>(£/) | F^Ei) F £f=£i ^ 0+ea (E f ) d£i 

Ei + ep — e a 



A = E f 



T(A,Ti,T f ) BW 

F{A,(Ji,(J f)Gauss = 



(20) 



2tt a 2 + (r i + r / ) 2 /4 

1 A 2 



exp 



2(af + aj) 



Eq. (20) with Tew was derived earlier j|, |i~9"| . The procedure used in deriving diag term is equiv- 
alent to starting with H = h(l) + V(2), first generating the non-interacting paricle (NIP) bivariate 
strength density I^ iv .Q(x, x') due to h(l) (using the results of (T^j), then convoluting the NIP spikes 
with independent spreading functions, due to V(2), at the two energies and finally simplifying as 
in (20). This assumes that the spreadings are uncorrelated. However, the of f diag term in (16) 
arises due to correlations and the bivariate correlation coefficient £ = (o^VOV^j / (O^Oj (VV), 
takes them into account, i.e. the spreading function Pb iv .Q{y,y') in 



l biv;0 



®pl v . [E,E') 



is not Pg.i(y)pg. f (y') but it is pX iv . (y, y'; 0, 0, o-i, o-f, C). Following JO 
with C the final result is same as in (20) but with a modified T, 



it is easily seen that even 



T(A 



A 2 



2vr(a 2 + 0- 2 - 2^07) 



exp 



2(af + a* 



f 



(21) 



Therefore for A = and £ — > 1, there is enhancement in the matrix elements compared to the diag 
approximation. This is indeed seen in EE(l+2) calculations earlier |l{| and now it is clear that 
with C in p\i v -o we have the proper theory for transition matrix elements of one-body operators 
(the theory with (21) is expected to operate in the chaotic domain defined by A > \p k and Xfu is 
defined in Section 5.1) . In fact for <j, ~ cry, the enhancement is — £). Note that for £ = 

we get back the diag approximation which is a GOE result. Also as noted in Jl^], the enhancement 
grows with m as £ grows with m (see (9)). Applications of the results in (20,21) and their variants 
[15], for nuclei and atoms are given in Jl5|, 20] and [0, 21] respectively. 



6. Conclusions 
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In this article an attempt is made to give an overview of the subject of embedded random matrix 
ensembles (classical EGOE and its various deformations) for complexity and chaos in interacting 
particle systems. It should be clear that generic embedded ensembles results are relevant in the 
quantum chaotic domain of isolated finite interacting many particle systems such as nuclei, atoms, 
molecules, atomic clusters, quantum dots etc. and therefore large scale explorations of deformed 
embedded ensembles are called for. 
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